Charge Density Wave Formation in the Low- Temperature-Tetragonal Phase of 

Cuprates 
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We calculate the influence of transverse fluctuations on the longitudinal dynamics in the striped 
phase of cuprates by using the bosonization technique. We find that a charge density wave instability 
can arise if the stripe is quarter filled and the underlying lattice potential has a zigzag symmetry. Our 
results explain why static stripes are experimentally observed in underdoped La2-i- !/ Nd H Sr a: Cu04 
exactly at the onset of the low-temperature-tetragonal transition. 
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During the last years, a large number of experiments 
have confirmed the formation of ID charge stripes in 
doped Mott insulators Gr9. Recently, it was shown 
H that for La 2 - 2; Sr 2 ;Cu04 co-doped with Nd the stripe 
structure becomes static and incommensurate peaks cor- 
responding to charge- and spin-stripe order can be mea- 
sured by quasi-elastic neutron scattering. Moreover, it 
was observed that for underdoped compounds (x = 0.10 
and x = 0.12) the temperature T co at which charge order 
sets in coincides with the onset of a structural transition 
from a low-temperature-orthorhombic (LTO) to a low- 
temperature-tetragonal (LTT) phase, TltTj whereas for 
the optimally doped compound T co < Tltt @- There- 
fore, although the formation of static stripes seems to 
be undeniably connected to the structural transition, it 
is not yet clear how and why this should happen, and 
moreover, why there is no correlation between T co and 
Tltt at optimal doping. 

In this letter, we address these questions and show that 
a charge density wave (CDW) instability can arise if two 
conditions apply simultaneously, namely, if the stripe is 
quarter filled and if the underlying lattice potential has 
a zigzag symmetry. This is the situation for underdoped 
cuprates in the LTT phase: neutron scattering experi- 
ments have shown that for (0.06 < x < 0.12) the stripes 
are incompressible, with one hole every two lattice sites 
HQ. In addition, within the LTT phase the buckling of 
the oxygen octahedra has the zigzag symmetry needed to 
stabilize bond centered stripes. The consequence of the 
combined presence of both facts leads to the formation 
of static stripes. 

We calculate the influence of transverse fluctuations on 
the longitudinal dynamics along the stripes. The longi- 
tudinal modes are described as Luttinger liquids. The 
coupling between the longitudinal and transverse modes 
arises because the electrons moving along the stripe are 
scattered from transverse kinks. By writing the Hamilto- 
nian in a bosonic representation, we find several forward 
scattering channels and one 4kp scattering momentum 
transfer from longitudinal to transverse modes. The lat- 



ter becomes relevant at quarter filling and leads to the 
formation of a longitudinal CDW if the transverse modes 
are frozen with a zigzag symmetry. 

Before introducing our model of a Luttinger liquid cou- 
pled to transverse kinks, we briefly review the transverse 
dynamics of a chain of holes. This problem was con- 
sidered by Zaanen et al. ||, who showed its equivalence 
to a quantum spin-1 chain. The full phase diagram for 
this problem was determined numerically by den Nijs 
and Rommelse @] after earlier calculations by Schulz [0 , 
who treated the spin-1 problem as two coupled spin-1/2 
chains. Here, we review Schulz calculations B and gen- 
erate the full phase diagram from this formalism. Then, 
we study the effect of a coupling between the longitudinal 
and the transverse sector to determine when the operator 
responsible for the CDW instability becomes relevant. 

Our starting point for the transverse modes is the spin- 
1 Hamiltonian |J 
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By replacing the spin-1 operator by a sum of two 
spin 1/2 operators, S n = T a (n) + T&(n), performing a 
Jordan- Wigner transformation into a system of interact- 
ing fermions, and then using a standard bosonization 
technique, one finally obtains H = H + + H- with 
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Here, *± = (*„ ± * 6 )/v/2 and x ± = {Xa ± Xb)/V2. The 
bosonic phase fields are 
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Xa/M = W [p R , a/b - p L , a/b ] e- a ^ 2 - ikx , (2) 

where a is the lattice cutoff and I the system 
size. The fields with different a/ b indices commute 
and [^ a /b{x),Xa/b{y)] = i5(x - y). The ^ a/b are 
related to density fluctuations of the a,b fcrmions 
through d x f a/b (x) = —K[p R<a / b (x) + p L , a/b (x) - 1/2]. 
The velocities of the acoustic excitations are u + ~ 
ts/i + 2(d-ZJ z )/{'Kt) and u_ ~ ty/1 - 2(d + J z )/(vrf). 
Further, K+ = [1 + 2(d - 3 J z ) / {nt)]- 1 / 2 , K_ = [1 - 
2(d + J z )/(nt)]- 1 / 2 and 9_(x) = n f dx'x-(x')- The 
coupling constants for the nonlinear terms are approxi- 
mately given by p\ = p% = J z + d and /13 = — t. While 
the p\ and pi terms represent, respectively, back- and 
forward scattering contributions, the p^ term originates 
from a non-local Fermi coupling and is the so-called disor- 
der operator. The most relevant operator resulting from 
the original Hamiltonian is cos (\/20_). There are other 
operators which mix the (+) and (— ) sectors but their 
scaling dimension is such that they are never dominant. 

To analyze the possible phases of the model defined 
through Eq. ([l]), we treat the non- linear terms as per- 
turbations and consider their scaling behavior Q. The 
Pi term is relevant for d > 3J Z . The transition at 
d = 3Jj is the roughening transition of the string, 
which can be seen calculating the large n behavior of 
the correlator of the transverse displacements G u (n) = 
(u(O)u(n)) = E"jio <SfSJ>> where S? = u i+l - Ui . 
Using the bosonized form S z (x) ~ — (\/2/Tr)d x ^+ + 
(2/Tra)e l7TX cos(-\/2*+) cos(\/2#_), we find the equal 
time correlation function G z = (S z (0)S z (x)) ~ 
{2K + /n 2 x 2 ) + C z e l7TX \x\~ 2K+ ~ 2K ~ which consists of a 
smooth and an oscillating part (C z is a constant). The 
smooth part gives rise to a logarithmic divergence of 
G u (n) for K + ^ 0. However, for d > 3J Z the p\ operator 
is relevant and K + scales to zero, killing the logarithmic 
divergence of G u (n). 

The (— ) sector is more involved. The P2 operator is 
relevant for d + J z < and /13 is relevant for d + J z > 
— 15irt/2 Thus, the gaussian fixed point is never sta- 
ble and the (— ) sector always flows to strong coupling. 
This implies that a perturbative approach starting from 
the Gaussian model cannot be trusted. However, it was 
pointed out by den Nijs [^J^] that the (— ) sector has an 
Ising symmetry. It has a form identical to the continuum 
limit of a classical 2D XY model H xy = Jt2uj) cos {4'i~ ( t ) j) 
with a 2-fold (Ising) symmetry breaking term of the form 
cos 2<pi. We therefore expect a transition in the (— ) sector 
of the Ising type. We locate the critical line of the tran- 
sition as the line where the scaling dimensions of the two 
non-linear operators in H- are equal, d + J z ~ —3tn/2. 

To analyze the properties of the phases we investi- 
gate the behavior of the equal time correlation function 
G ± (n) = (S+(n)S-(0)). If *+ is ordered (i.e. d > 3J Z ), 



then G± decays exponentially as 8+ correlations are 
short ranged. Similarly, if d < 3J Z and d+ J z > —3tn/2, 
such that has long range order, then 0_ correlations 
also are short ranged, again leading to an exponentially 
G±. However, for d < 3 J z and d+J z < —itn/2 : 0_ is or- 
dered and G± decays algebraically, G±(x) ~ |a;| _1 / 4 ^+. 
Combining the results from the (+) and (— ) sectors, we 
then find five different phases. 1) a gap- full flat phase, 
with dominant pi and /13 perturbations, exponentially 
decaying G z and G± correlations and a finite limit of 
G u . 2) A gap- less rough phase in which only p 3 is rele- 
vant, with an algebraic decay in the smooth part of G z , 
exponentially decaying G± and logarithmically divergent 
G u correlations. 3) A gap-less bond centered (BC) rough 
phase with a zig-zag pattern (anti-ferromagnetic correla- 
tions with S z = ±1 but no S z — states). This phase dif- 
fers from the rough phase in that it has an algebraic decay 
in the G± correlations. 4) A BC flat phase with domi- 
nant pi and P3 perturbations. Like the BC rough phase, 
this phase has a zig-zag pattern which now is however 
long ranged. This phase is gaped. 5) For 3 J z — d> wt/2, 
the present analysis cannot be applied because the di- 
agonalization of the quartic a, b interactions through a 
Bogoliubov transformation breaks down, leading to an 
unphysical purely imaginary value of K + . As argued by 
Schulz this signals a transition to a ferromagnetic 
state, i.e. a diagonal stripe state. 

Besides the five phases found by Schulz, den Nijs 
has identified a sixth phase, see Fig. 1. The additional 
phase found by den Nijs is a disordered flat phase (DOF), 
which is gap- full. This phase has, in contrast to the flat 
phase, a finite density of kinks and anti-kinks (S z = ±1 
states) but the G u correlator does not show a logarithmic 
divergence, making this phase different from the rough 
phase. In the DOF phase the kinks are positionally dis- 
ordered, but have an anti-ferromagnet order in the sense 
that a kink S z = 1 is more likely to be followed by an 
anti-kink S z = — 1 rather than another kink, with any 
number of S z — states in between them. In spin lan- 
guage, the DOF phase is the valence bond phase which 
is responsible for the Haldane gap. 

In fact, there are signs of this transition also in the 
Abelian bosonization approach. At d + J z = 0, p\ and 
pi changes sign. This is unimportant for pi as in this 
parameter regime the (— ) sector is dominated by the /13 
and not the pi operator. However, if the p\ operator 
is relevant, then a sign change of p\ is of consequence. 
Typically, such a sign change is indicative of a competi- 
tion between two fixed points. A similar situation occurs 
in a spin- 1/2 chain with nearest and next-nearest inter- 
actions, where dimcrized and anti-ferromagnetic ground 
states compete |l0| ]. In the present situation, the com- 
petition is between on-site (d) and nearest neighbor in- 
teraction (J z ). A useful order parameter to distinguish 
between the flat and the disorder flat phases is the parity 
of the steps, P — (cos(7rii„)), which in bosonized form 
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becomes P ~ (cos [-y/2 (*+(x) - *+(0))]), with x = 
being the boundary of the string. If the p\ operator is 
relevant, the field i&+ is pinned and a sign change in ji\ 
leads to a 7r phase shift of v / 8*+- As a result, P shifts 
from zero to a finite value (if the boundary term is left 
unchanged). Therefore, while the Abelian bosonization 
approach does not allow for a detailed study of the pre- 
roughening transition or the DOF phase, the existence of 
this transition can be readily inferred. 
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FIG. 1. Phase diagram for the spin-1 chain and the cor- 
responding stripe configurations. The inset in the left corner 
shows a 4 }cf CDW state in the BC flat phase. The size of the 
circle represents the local density of holes, a) commensurate 
ground state, b) CDW domain wall soliton. c) soliton state 
of both the transverse and longitudinal order. 

We now study a model including both, transverse and 
longitudinal modes. The longitudinal modes along the 
stripe are described by a Luttinger liquid |1 1|, 
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where g\ is the backward scattering strength. The fields 
<f p/r7 = ($ T ±$;)/V2, U p/a = (n T ±IIj)/V2 are related 
to the Fermi fields of left- and right-moving particles. 
The parameters of the bosonic model depend on the un- 
derlying lattice model. We assume repulsive interactions, 
therefore K p < 1 and K a > 1. The backscattering term 
is irrelevant under the renormalization group procedure 
and will be neglected here pT| . 



Physically, the coupling of the transverse and longitu- 
dinal modes arises from the scattering of the electrons 
which move along the stripe from transverse kinks. This 
interaction is local and we can express it in terms of the 
longitudinal electron density p^ n and the kink density 
S 2 . We will consider only terms that couple the charge 
density to the kink density. Terms that couple the kinks 
and the spin density are irrelevant in presence of repul- 
sive interactions and we thus omit them. A local coupling 
of the electron density to the kink density can thus be 
written as 
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with coupling constants 71,72- These terms can be read- 
ily bosonized. From H c i one obtains the following terms, 
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with Ai = 7i/ (\/2ir 3 a 2 ). The first term represents for- 
ward scattering between the acoustic transverse and lon- 
gitudinal charge modes. In general, these two modes have 
different velocities u + 7^ u p and thus this interaction is 
retarded as can be seen by performing a Gaussian av- 
erage over one of the modes. This first term and the 
Gaussian parts of H+ and H p and can be jointly di- 
agonalized, leading to a hybridization of the transverse 
and longitudinal modes. This leads to an instability of 
the system which is similar to the Wentzel-Bardeen in- 
stability in one dimensional metals coupled to phonons. 
However, this instability occurs only at a very large cou- 
pling 71 - 4y/(nK + K p ) I {u p u+) @. For small 71, the 
case considered here, the zero momentum transfer inter- 
action is not very efficient |D||. Also, as the hybridization 
is small for small 71, the scaling analysis below is only 
weakly affected by the hybridization and we thus ignore 
this correction. The Ai and A2 terms are only impor- 
tant if the respective cosine terms have a finite expec- 
tation value. If they do, these terms act like a shift of 
the chemical potential of the electrons, as d x Q p measures 
the deviation from the charge density of its equilibrium 
value. These terms compete however with terms gener- 
ated by 72, as will be discussed below. Bosonization of 
the 72 interaction gives the following terms, 

H C 2 ~ j dx A 3 d x ^+ cos (y8$ a ^ + A 4 cos (y8$ a ^ x 
cos (%/81' + ) + A 5 cos (y&Qp - a/2*+ + (4fc F - 7r) x\ x 

(6) 
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The A3 and A4 terms result from scattering involving 
two a,b operators whereas A5 involves four a,b operators. 
While the first two terms describe forward scattering, the 
A5 term results from processes with momentum transfer 
4fcp from longitudinal to transverse modes. At kp = 7r/4, 
i.e. at quarter filling of the stripe this term is important, 
as the oscillatory x dependence vanishes. 

To understand the effect of the longitudinal-transverse 
coupling we examine the scaling dimensions of the Ai 
operators. We envision strongly repulsive interactions 
among the electrons on the stripe with K p < 1 and K a > 
1 . The A4 operator is relevant only for K„ + K+ < 1 and 
thus unimportant. Similarly, A3 can be neglected as it is 
relevant only for K a < 1/2. The Ai operator is relevant 
for A"_ < 1/2, A2 is relevant for K + < 1/2 and the A5 
operator is relevant for 2K p + K + + A_ < 2. 

Let us consider that K + < 1, i.e. the (+) sector is 
massive. In that case, the transverse modes behave at 
large wavelengths as if K + = 0. Then, if 2K p + A_ < 2, 
the A5 operator will become relevant at large scales. The 
A5 term pins the longitudinal charge modes to the trans- 
verse modes and induces a longitudinal CDW. This is 
seen for the case with A'_ < 1/2, i.e. in the BC flat 
phase, where both ty + and have long range order. 
The freezing of the transverse modes in a zig-zag pat- 
tern leads to a t periodic potential for the transverse 
modes and allows for Umklapp scattering of the quar- 
ter filled stripe which causes a longitudinal 4fcp CDW. 
Thus, a relevant A5 term makes the $ p field massive and 
pins it to the frozen kink-antiking pattern of the trans- 
verse fluctuations. Such a stripe with a frozen kink pat- 
tern and a longitudinal CDW pinned to it is depicted in 
Fig. 0a. If the A5 term is relevant, it is in direct compe- 
tition with the Ai ; 2 terms which favor a non-zero d x & p . 
While an exact treatment of this competition is not cur- 
rently possible, we can gain some insight if we assume 
mean field values for the transverse fields. Deep into 
the BC flat phase we may replace cos (\/8^±) in the 
Ai,2 terms by (cos (V& i &±))- Similarly, the A5 term be- 
comes A 5 / dx cos (\/8$ p ) (cos (V2*+)) (cos (V2*-)). 
In this approximation the longitudinal charge sector 
has a Hamiltonian which is well known from studies of 
commensurate-incommensurate transitions [|l4[ , p~5| . Such 
a model can be studied within a fermion description in 
terms of spin- less holons 0,0), in which the A5 Umk- 
lapp term becomes quadratic rather than quartic as in 
the original Fermi operators. In these models, a gaped 
commensurate state, which exists for small A^a (i.e. a 
small prefactor of the d x ^ p term) is followed by an in- 
commensurate state at large values of Ai 2 where the gap 
is destroyed and a small number of holons move as weakly 
interacting particles. These holons are solitonic in char- 
acter and can be interpreted within the present frame- 
work as point like domain "walls" of the it periodic CDW. 
Such a CDW soliton state, which involves only the trans- 



verse sector, is depicted in Fig. [l|b. A different soliton is 
also possible, see Fig. |l|c, which represents a domain wall 
for both, transverse and longitudinal degrees of freedom. 

It is interesting, that even if A_ > 1/2, i.e. in the 
DOF phase, the A5 term can be relevant for sufficiently 
strong repulsion (small A p ). However, if one assumes 
a Hubbard model for the longitudinal modes, the lower 
limit for K p is 1/2 so that in the flat phase, with A_ > 1, 
the A5 operator is always irrelevant. This result has two 
important consequences: first, since A5 is irrelevant in 
the flat phase, any kind of underlying structure which 
favors site-centered stripes is unable to produce a CDW, 
implying that the usual stripe "cartoon" Q is actually 
misleading. Second, it allows us to understand why a 
static CDW is formed within the LTT phase of under- 
doped cuprates, since the symmetry of this phase favors 
a bond-centered pattern, for which the A5 operator is 
relevant. 

In conclusion, we evaluated the effect of transverse fluc- 
tuations on the longitudinal transport in a Luttinger liq- 
uid stripe phase. We found that a CDW instability arises 
if the stripes are quarter filled and the underlying po- 
tential has a zigzag symmetry. Our results explain why 
static CDWs are formed in under-doped cuprates exactly 
at the onset of the LTT structural transition. 

C. M. S. was supported by Swiss National Foundation 
under grant 620-62868.00. 
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